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ABSTRACT 

Current microlensing searches calibrate the mass fraction of the Milky Way halo 
which is in the form of Massive Compact Halo Objects (MACHOs). We show that 
surveys like the Sloan Digital Sky Survey (SDSS) can probe the same quantity in halos 
of distant galaxies. Microlensing of background quasars by MACHOs in intervening 
galaxies would distort the equivalent width distribution of the quasar emission lines by 
an amplitude that depends on the projected quasar-galaxy separation. For a statistical 
sample of ~ 10^ quasars as expected in the SDSS, this distortion is detectable at the 
2 cj level out to a quasar-galaxy impact parameter of several tens of kpc, as long as 
extragalactic halos are made of MACHOs. Detection of this signal would test whether 
the MACHO fraction inferred for the Milky-Way halo is typical of other galaxies. 


Subject headings: cosmology: dark matter — galaxies: halos 
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1. Introduction 

One of the basic unknowns in cosmology is the fraction of the dark matter mass which is 
in the form of compact objects, such as stellar remnants. Paczyhski (1986) proposed monitoring 
of background stars in the Large Magellanic Cloud (LMC) and searching for rare amplification 
events due to microlensing as means of constraining the mass fraction of compact objects in the 
Milky Way halo out to 50 kpc. Current microlensing surveys suggest that ~ 60 ^( 20 % of the mass 
of the halo of the Milky Way galaxy might be in the form of Massive Compact Halo Objects 
(MACHOs) with a mass Mq (Alcock et al. 1996). It remains to be seen whether future 

studies confirm these early reports. 

The search for brightening of individual background stars due to microlensing by foreground 
MACHOs cannot be trivially extended beyond the scale of the Milky Way halo. Images of galaxies 
at large distances include many faint stars within each resolution element (pixel), most of which 
are unlensed. Alternatively, one can search for a weak temporal fluctuation due to the chance 
amplification of one of the stars within a given pixel (Gould 1996; Crotts 1996). This method is 
currently being applied to the nearest massive galaxy, M31 (Crotts & Tomaney 1996). 
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At yet larger distances, MACHOs can be detected through their effect on quasars. The lensing 
cross-section of a MACHO increases in proportion to its distance from the observer. Thus, despite 
the fact that quasars are rare, the likelihood of a quasar which happens to be located behind an 
extragalactic halo to be microlensed is ~ 5-6 orders of magnitude larger than that of an LMC star 
behind the Galactic halo. This makes quasars useful sources for the study of extragalactic halos 
as LMC stars are for the Galactic halo. The only question is how can one identify the existence of 
microlensing in the former case. 

The crossing time of the lensing zone of solar mass MACHOs in the Galactic halo is of the 
order of a month, and so microlensing of LMC stars can be conveniently identified through the 
achromatic variability signal of an otherwise steady stellar source. However, if the same halo 
objects are placed at cosmological distances, the duration of their lensing event is stretched by 
~ 2-3 orders of magnitude up to the inconvenient scale of decades. Moreover, the lensing events 
of quasars would be contaminated by noise from the uncertain level of intrinsic variability of 
the quasars on the same timescale (Maoz 1996). Therefore, variability is not an efficient search 
technique for microlensing by solar mass MACHOs in distant halos.|^ 

An alternative method that does not require monitoring over time is to search for a change in 
the equivalent width of the broad emission lines of the quasar. The characteristic Einstein radius 
of a solar mass lens at a cosmological distance is ~ 5 x 10^® cm, comfortably in between the scales 
of the continuum-emitting accretion disk (;^ 10^® cm) and the broad line region (~ 3 X 10^^ cm) of 
a bright quasar. This implies that a microlensing event would significantly amplify the continuum 
but not the broad lines emitted by the quasar. As a result, the equivalent width distribution of 
the broad emission lines (Francis 1992) will be systematically distorted in a sample of microlensed 
quasars (Canizares 1982, 1984). Dalcanton et al. (1994) used the lack of redshift evolution in the 
equivalent width distribution of quasars to limit the mean density parameter of MACHOs with 
a mass 10“^-60 Mq in the Universe, ^ 0.2. This limit, however, does not exclude the 

possibility that MACHOs account for most of the dark matter in galactic halos (or even globally, 
if the Universe is open). Since luminous stars cluster in galaxies, it is possible that their dim 
counterparts are also concentrated around galaxies. 

The probability for microlensing depends on the projected separation between the background 
quasar and the center of the intervening galaxy. In this paper we propose to look for MACHOs 
in extragalactic halos, by searching for the equivalent width distortion of quasar emission lines as 
a function of the projected quasar-galaxy separation. The forthcoming Sloan Digital Sky Survey 
(SDSS; see jhttp:/ /www.astro.princeton.edu/BBOOK|) , will catalog ~ 10^ quasars and ~ 10® 
galaxies, and would provide an ideal data base for such a study. 

The outline of the paper is as follows. In §2 we describe our model assumptions. In §3 we 


^Although Jupiter mass lenses could cause quasar variability on the more convenient timescale of <year (Hawkins 
& Taylor 1997), current microlensing searches (Alcock et al. 1996, Ansari et al. 1996) rule out the possibility that 
such objects contribute significantly to the mass of the Galactic halo. 
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derive the expected distortion of the equivalent width distribution of background quasars as a 
function of the MACHO mass fraction in the intervening galactic halos. In §4 we quantify the 
expected lensing signal in SDSS, using Monte-Carlo realizations of the quasar and galaxy fields. 
Finally, §5 summarizes our main conclusions. 


2. Model Assumptions 


We model a galactic halo as a singular isothermal sphere (SIS), whose average surface 
mass density at a distance ^ from the center is S(^) = ct^/( 2G'^). The one-dimensional velocity 
dispersion, cr, is related to the luminosity, L, of the galaxy through the Tully-Fisher (1977) 
relation, u/it* = (L/L*)", with a ~ 0.4 in the R band (Strauss & Willick 1995), cr* ~ 170 km s“^ 
and L* = 10^*^ LQh~‘^. It is convenient to normalize all length scales in the lens plane by 


Co = 47r 


a 




( 1 ) 


where Doi, Dos and Dig are the observer-lens, the observer-source, and the lens-source angular 
diameter distances, respectively. These distances depend on the cosmological parameters; in the 
following we assume H = 1, A = 0 and h = 0.5. Dehning x = C/Co as the normalized image 
position in the lens plane, the amplification factor of a point source by a SIS is 

X 


(/^) = 


a: — 1| 


( 2 ) 


Multiple images occur when x < 2. For x > 2, the single image position is at x = y -|- 1, where 
y = f/Doi/(DosCo)) and rj is the distance of the source from the optical axis in the source plane 
(Schneider et al. 1992). 


Equation (|^ yields the amplification factor ji only in the case where the mass distribution is 
perfectly smooth. However, in the presence of MACHOs one obtains a probability distribution 
of amplification factors at any given value of x. The average amplification with this distribution 
must still be given by equation (^) to guarantee flux conservation. Since the lensing zone of a 
solar mass MACHO is six orders of magnitude smaller than the length scale over which the global 
S varies, the MACHOs can be treated locally as a field of randomly distributed point masses with 
some average surface density S. We then define p{iJ., x, e)dfi as the probability that a point source 
whose line-of-sight passes at an impact parameter x relative to the center of a galaxy is amplified 
by a factor between fj, and p + dy. This probability depends on the MACHO mass fraction e, which 
we take as a constant. Schneider (1987) and Bartelmann &: Schneider (1990) derived the following 
analytical expression which agrees reasonably well with results from numerical simulations, 


p{y,x,e) 


'C X {p - /i*)^ 
. /D+3 


-|- A X (/r 


/i*) exp(—i3/i)d/i 


H{p — p^)dp. 


( 3 ) 


Here denotes the Heavyside step function; A, H, and C are functions of x -the first two 

are determined by the constraints of normalization and flux conservation and the third by the 
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limit of high amplification (Schneider 1987). The parameter r/ is chosen so as to make A and 
B well-defined and positive for all possible values of x. The value of r] rises monotonically with 
decreasing e, because the probability distribution approaches a (5-function as the MACHO fraction 
goes to zero. The parameter is a function of e such that it approaches (/u) in the limit of e ^ 0 
(see Bartelmann &: Schneider 1990 for full details). 

Our goal is to quantify the distortion of the equivalent width distribution of quasar emission 
lines due to microlensing by MACHOs in halos of foreground galaxies. For a given magnitude- 
limited sample of galaxies and quasars, we need to select the galaxy and quasar properties at 
random. We assume that the local luminosity function of galaxies has the Schechter form, 

<t.a(L)dL = 4'. «?{-£}<'(!;). W 

with = 0.018 {h/0.5)^ Mpc~^ and zz = 0.3 in the R band (Lin et al. 1996). Since a typical 
survey contains galaxies at very low redshifts (e.g. (z) ~ 0.1 for the SDSS), we ignore evolutionary 
corrections and /c-corrections to this luminosity function and scale the density of galaxies in 
each luminosity bin by (1 -|- z)^. On the other hand, the quasar density peaks at high redshifts 
z = 2-4 (see Fig. 2 in Shaver et al. 1996). Because of the disparity between the lens and source 
redshifts, we can assume that all the quasars are at an effective redshift Zg ~ 3 and make use 
of their well-determined total number count. The inaccuracy introduced by this approximation 
to the lensing calculation is small (~ 10-20%) and not significant relative to the observational 
uncertainties in the redshift evolution of the quasar luminosity function. In this approach, we 
only need to use the integral number counts of quasars as a function of B band flux, which is well 
described by a broken power-law 


N{> S) = C' X 


(S/So)-^^ if5<5o 

(5/5o)“^" if5>5o 


( 5 ) 


where (3i ^ 1 and P 2 ~ 2.5 (e.g., Setti, 1984; Marshal, 1985), and C ^ 1/deg^ (Hartwick & 
Schade 1990). The break flux Sq corresponds to a B magnitude ~ 19.5. We can safely ignore any 
correlation between the quasars and the galaxies because of the clear separation between their 
redshifts. 


3. Signature of Microlensing 

The flux of a quasar can be significantly amplified by microlensing only if the size of its 
emission region is smaller than the projected Einstein radius of the lens in the source plane, te. 
The maximum amplification of a circular source of radius Vg and uniform brightness is given by 
[Schneider et al. (1992), p. 38] 


Mmax 


Y^lT4(rE/r^ . 


( 6 ) 



- 5 - 


At cosmological distances, the Einstein radius of a star of mass Mgtar obtains the characteristic 
value te ~ 5 X 10^® (Mstar/Af©)^/^ cm. In comparison, the optical continuum emission of quasars 
is believed to originate from a compact accretion disk. The UV bump observed in quasar spectra 
is often interpreted as thermal emission from an accretion disk with a surface temperature 
Tdisk = 10^ Ts K, where Ts 1 (e.g. Laor 1990), and so the scale of the disk emission region 

must be ~ 10^^ cm x , where L^q is the corresponding luminosity of the quasar in units 

of 10^® erg s“^. Thus, for lens masses Mgtar S> 10~^Mq, the continuum source is much smaller 
than the Einstein radius and could therefore be amplified considerably. This expectation is indeed 
confirmed in the nearby lens of Q2237+0305, where variability due to microlensing has been 
observed (Wambsganss et al. 1990; Rauch &: Blandford 1991; Racine 1991; see also Gould &: 
Miralda-Escude 1996). On the other hand, reverberation studies of the time lag between variations 
in the continuum and the line emission in active galactic nuclei indicate that the broad emission 
lines of quasars originate at a distance of ~ 3 X 10^^ cm (Peterson 1993; Maoz 1996; Netzer 
&: Peterson 1997). Eor a solar mass lens and a quasar with L^q ~ 1, equation (|^ implies that the 
maximum line amplification differs from unity by only ^ 10%. This argues that microlensing of 
the broad line region by a single star can be neglected. The broad line region would, however, 
be macrolensed by a factor (fi) due to the average effect of the galaxy as a whole. As a result of 
this differential amplification of lines and continuum, the equivalent width of the lines will change 
during a microlensing event. 


More specifically, let us define /j/(Ao) to be the intensity of the quasar continuum in the 
neighborhood of the wavelength Aq of a particular emission line and A/i^(A) to be the difference 
between the total (line+continuum) intensity and the continuum. If the continuum is amplified by 
a factor //, then its observed intensity changes to The intensity of the lines is enhanced by 

a factor (//), given by equation (^). Consequently, the equivalent width (EW) of the emission line. 


W = 


AI^ 


dX 


is changed by a factor (^)//r, namely 


W = wM 


where Wq is the intrinsic EW of the unlensed quasar. 


(7) 

( 8 ) 


Even in the absence of lensing, quasars do not possess a single EW value in their emission 
lines but rather show a wide probability distribution of EW values (Erancis 1992), which we define 
as P{Wo) and model by a log-Gaussian form. 


P{Wo) 



exp 


-[In(IEo) - 



(9) 


The parameters uj and aw obtain different values for different emission lines. The distribution in 
equation i) is assumed to be observed in the absence of lensing. However, if the line of sight to a 
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quasar passes at an impact parameter x from the center of a galaxy with a MACHO fraction e, 
then its equivalent width is drawn from the distribution 

P{W,x,e)=J^ dnp{n,x,e)P (10) 


Figure 1 shows the probability distribution in equation (1C) as a function of x, for e = 0.3 
(panel a), e = 0.5 (panel b) and e = 1 (panel c). The unlensed distribution was selected for the 
Mgll line, based on the data in Francis (1992). The microlensing signal is similar for other lines 
(Perna &: Loeb 1997). Naturally, the systematic distortion of the lensed distribution increases as 
X decreases, because of the corresponding increase in the surface density of MACHOs. However, 
the number of observed quasars declines as x becomes smaller, and so the signal suffers from an 
increasing statistical noise in this limit. In order to properly evaluate the signal-to-noise ratio we 
need to combine the above effects. 


Let us consider a galaxy with luminosity L at redshift z and a population of background 
quasars with a flux > 5 qso and an average number per solid angle N{> S'qso)- We define 6s to 
be the angle between the line-of-sight to a given quasar and the optical axis of the galaxy. The 
number of quasars within a ring of width dOg around 6s is N(> SQso)2vr0^d0s- Using 6s = rj/Dos, 
and defining 6^ = ^o/Dd, this number can be expressed as 

N{> 5Qso)27r^ = N{> SQso)‘27r^ydy = iV(> 5qso, x, e)27r02(x - l)dx. (11) 

^OS ^o\ 


The last step used the relation between x and y for a single image in the SIS model and 
incorporated amplification bias through the dependence of the local density of quasars with 
observed flux > Sqso on x, 

N{> SQso,x,e) = j .. f dpp{p,x,e)N{> (12) 

(/i(x)) Ji p 

Equation (|^) gives the number of quasars per unit solid angle inside a ring of radius x and width 
dx around the center of a single galaxy. To obtain the total number of quasars per solid angle 
which are seen within a distance dx of x of all galaxies up to a limiting flux S'cah one has to 
integrate over the galaxy distribution, both in luminosity and in redshift. 


N{> S'qso, > Seal, x, e)dx = 


X 


27r{x — 

l)dx 

(pix)) 

rzs 

dt 

/ dz 

c— 

Jo 

dz 


^ dp p{p, x, e)N 

roo 

(l + z)^ dL 4 >g{L)^q{L, z). 


(13) 


Here we have implicitely assumed that the rings around different galaxies do not overlap. As 
shown later, the signal-to-noise ratio of the microlensing signature peaks at small values of x, 
where this assumption appears to be well satisfied. If Soai is the minimum flux needed to detect a 
galaxy, then Loai in equation (|l^) is given by Loai = 47rZl|(z)SGai, with the luminosity distance 
Di^{z) = (1 z)‘^Do\{z). 
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4. Applications and Resnlts 


Next, we apply the results from §3 to simulated data for a survey similar to SDSS. Let 5 qso 
and Scai be the minimum fluxes needed for the detection of quasars and galaxies, respectively, 
and let Afl be the total solid angle surveyed on the sky. By definition, the overall (sky averaged) 
source count of quasars is not influenced by the amplification bias. Hence, the total number of 
detected quasars is given by Atot = N{> Sqso)^^^ with N{> S) from eqnation (^). 

To simulate the results of real observations, we start by generating Atot random numbers 
drawn from the distribution 


P{x)^ 


A(> Sqso) 


A(> Sqsoj > ‘S'cab X, e) 


(14) 


such that 


dx N{> Sqso, > •S'cab x, e) = N{> Sqso) 


(15) 


with N{> Sqso,> SGa.hX,e) given by equation (|^). In onr actnal calculation, the value of the 
lower limit of integration has been set to 2 to avoid the complications arising from multiple images. 
Since the fraction of quasars which are multiply imaged is very small (^ 1%), our statistical 
analysis is not affected by this assumption. The value of Xmax is fixed by equation (IT) to reflect 
the crude boundary beyond which the halos of foreground galaxies overlap; its precise value has 
little impact on onr conclusions since the signal-to-noise ratio is high only for values of x <C Xmax- 


We have generated a source catalog at random based on equation (|l^ , with equivalent widths 
W chosen according to the probability distribution in equation { 0 )- This data set is chosen 
to represent a mock realization of SDSS, for which we evaluate the signal-to-noise ratio of the 
microlensing signature assuming a spectral resolution for the lines of ~ 5A. Our analysis of this 
catalog involves several steps. First we bin all sources based on their equivalent width, Wi, and 
derive the “average” probability histogram Pave(^i)) which reflects equation (^) to within the 
statistical noise. We then bin the sources based on their lowest value of x relative to a foreground 
galaxy. (Note that the association of a galaxy to a quasar based on the minimum x-value does 
not necessarily pick the closest galaxy to the quasar on the sky.) The equivalent widths of the 
quasars in each Xj-bin are then binned, and a analysis is performed to test the significance of 
the deviation of the data in each bin from the average distribution Pave(kFi). More specihcally, 
we dehne ~ /rii], where hi is the number of quasars in the equivalent-width bin 

i within a given x-bin and n* is the expected number according to the average distribution Pave- 
The end point of each x^-bin is varied (starting from the end of the previous Xj_i bin) so as to 
minimize P{x^)- The smaller P(x^), the more significant is the deviation of the data in each bin 
from Pave- 


Fignres (2a), (b) and (c) show P(x^) as a function of the binned x for simulated sets of data 
with Zs = 3, and magnitude detection thresholds of m-Gai ~ 22 P-mag and ruQso ~ 22 P-mag, 
giving a total quasar count of Ntot ~ 10^, as expected for the forthcoming SDSS. In panel (a) the 




MACHO fraction is e = 0.3, in panel (b) e = 0.5, and in (c) e = 1. The first case corresponds to 
the situation where the halo dark matter is dominated by a smooth dark matter component (e.g., 
elementary particles), and the last case to a situation where all the dark matter is in the form 
of MACHOs. The difference between the three panels is apparent. For a value of e < 0.5, P{x‘^) 
is of order of several tenths, implying that the microlensing signal is not statistically significant. 
However, for e ~ 0.5, the value of P{x‘^) small x indicates that the EW distribution starts 
to become statistically different from Have as a result of microlensing. The difference between 
these distributions is strong when the MACHO fraction is close to unity. Figure (2c) shows that 
the equivalent widths of quasars whose lines of sight intercepts a galaxy at an impact parameter 
re ^ 10 from its center, can be rejected as being drawn from the distribution Pave with a confidence 
level 2cj (the threshold level indicated by the dashed horizontal line). We have found this result 
to be robust using several realizations of the mock SDSS catalog. 

The length scale associated with a given value of x = depends on the luminosity L and 
redshift zi of the lensing galaxy, and on the source redshift Zs', for Zs = 3, 


^0 = 96 


(1 + Z\) 


3 

.2-y/l + z\ 


1 

1 + zi 



(16) 


For example, an L* galaxy obtains = 3.7 kpc at zi = 0.5. Therefore, detection of a microlensing 
signal at X ~ 10 would probe galactic halos out to several tens of kpc, similar to the local 
microlensing searches. Since this scale extends well beyond the luminous cores of galaxies, our 
proposed method would sample the dark matter content of the intervening halos. 


5. Conclusions 

We have shown that the MACHO fraction in extragalactic halos can be inferred from the 
distortion in the equivalent width distribution of background quasars as a function of their 
separation from foreground galaxies. This method could be used to test whether the properties of 
the Milky Way halo as inferred from local microlensing surveys are characteristic of other galaxies. 

For the magnitude limit and survey area of SDSS, we find that if galactic halos are made of 
MACHOs, the ratio between the microlensing signal and the statistical rms noise is 2 out to a 
halo radius of several tens of kpc (see Fig. 2). Larger surveys might be able to probe the signature 
of halos with a smaller mass fraction of MACHOs. 


We thank Matthias Bartelmann for providing us with a computer routine which calculates 
the probability distribution of amplifications. This work was supported in part by the NASA ATP 
grant NAG5-3085 and the Harvard Milton fund (for AL) and by a fellowship from the university 
of Salerno, Italy (for RP). 
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Fig. 1.— The equivalent width distribution for the Mgll emission of quasars. The solid line reflects 
the assumed unlensed distribution. The other curves show the distributions expected for a sample 
of microlensed quasars whose line of sight passes at an impact parameter x = 2 (dotted line), 
X = A (dashed line), x = 6 (long-dashed line), and x = 8 (dotted-dashed line) from the center of a 
foreground galactic halo. In panel (a) the MACHO fraction is e = 0.3, in panel (b) e = 0.5, and in 
(c) e = 1. 
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Fig. 2.— The x^-probability as a function of the binned x for a simulated SDSS catalog in which 
the equivalent widths and the quasar-galaxy separation x of ~ 10^ quasars are measured. The value 
of P{x^) indicates the probability that the EW data in each x-bin are described by the average 
distribution Pave (IE*)) obtained by considering all the EW’s, with no reference to quasar-galaxy 
separations. In panel (a) e = 0.3, in panel (b) e = 0.5, and in (c) e = 1. 



